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SYSTEM FOR HIGH EFFICIENCY SPACECRAFT ORBIT TRANSFER 

BACKGROUND OF THE INVENTION 
[0001] The present invention relates to orbit transfer techniques, and more particularly to 
5 systems and methods for efficiently transferring a spacecraft from one orbit to another orbit. 

[0002] Typically, a spacecraft is launched into an initial orbit that differs from the final 
orbit that the spacecraft is to follow during its lifetime. A spacecraft launch vehicle places 
the spacecraft in its initial or launch vehicle insertion orbit. The launch vehicle, however, 
generally does not have sufficient capability to place the spacecraft into its mission orbit. As 
10 a result, the spacecraft must be transferred from its initial insertion orbit to its final mission 
orbit, typically by using a spacecraft on-board propulsion system. 

[0003] For a geosynchronous orbit (GEO) spacecraft, a velocity change (Delta- V) of 
approximately 1600 m/sec is needed to transfer the spacecraft from a typical injection orbit to 
a final GEO orbit. As one skilled in the art will appreciate, the time allocated to transfer the 
15 spacecraft to the final mission orbit is limited, so the spacecraft typically requires a high- 
thmst propulsion system to generate the needed velocity change for a large spacecraft dry 
mass (e.g., 2000 to 3500 kg). 

[0004] The high-thrust propulsion system can comprise a liquid apogee engine (LAE) that 
bums a mixture of hydrazine and oxidizer to generate a thrust of about 100 to 150 Ibf with a 

20 mass specific impulse (ISP) of about 325 seconds (the ISP is the ratio of the thrust to the 

mass flow rate, a higher value indicating that less fuel is needed to achieve the same applied 
impulse). The orbit transfer then can be performed by firing a sequence of four to six LAE 
maneuvers at or near the orbit apogee, with each maneuver lasting less than 1 hour. Using 
this approach, the transfer to the GEO mission orbit can be accomplished in less than 2 

25 weeks. However, the performance of this approach, as measured by the mass delivered to 
orbit on a given launch vehicle, is limited by the capability of the launch vehicle and the 
amount of chemical propellant that can be stored on-board the spacecraft. 

[0005] The LAE orbit transfer approach provides reasonable efficiency and a rapid orbit 
transfer, but for heavy spacecraft the combined launch vehicle and on-board LAE system can 
30 be insufficient to transfer the spacecraft to the mission GEO orbit. To address this problem, a 
hybrid orbit transfer approach recently has been introduced that performs the orbit transfer is 



in two phases. The first phase includes some number of LAE firings to transfer the 
spacecraft fi-om the initial injection orbit to an intermediate orbit, which typically has its 
perigee above the Van Allen radiation belts to limit solar array degradation during the final 
orbit transfer phase (phase 2). 

5 [0006] During phase 2, high-efficiency low-thrust thrusters, such as ion thrusters, are fired 
for a period of several weeks to several months to complete the transfer to the mission GEO 
orbit. The thrusters can be gridded ion thrusters or Hall Current Thrusters (HCTs), which 
generate thrust by ionizing xenon atoms and accelerating them through a potential gradient. 
The ion thrusters have low thrust (e.g., less than 0.1 Ibf), but have high specific impulse of 
10 fi-om about 2000 to about 3500 seconds. Because this phase uses low thrust thrusters, the 
thrusters must be fired for extended time periods to provide a significant benefit. 

[0007] The time available to perform the final orbit transfer phase (phase 2) depends on the 
orbit time that is acceptable to the customer. Because the customer almost always wants 
revenue generating service to begin as soon as possible after the spacecraft launch, the 
15 allocated time usually is not more than 90 days. The longer the allocated time, the longer the 
ion thrusters can fire, and therefore, the larger the achievable payload mass that can be 
delivered to the mission orbit. 

[0008] According to prior- art implementations of the above described orbit transfer 
approach, the ion thrusters are fired continuously while the spacecraft attitude is controlled so 

20 that the thrust vector follows a pre-specified inertial trajectory, referred to as a thrust vector 
trajectory. This approach is referred to as a continuous firing strategy. The thrust vector 
trajectory can be determined by numerical optimization techniques well known to those 
skilled in the art. For example, the continuous firing thrust vector trajectory can be calculated 
to provide a transfer from an initial orbit to a final orbit in the minimum possible time. The 

25 resulting minimum-time transfer-orbit trajectory consists of a specific time history of the 
orbit elements as they vary fi-om an initial orbit to a final orbit while the thrusters are firing 
continuously. With this strategy, continuous thruster firing typically is maintained, except in 
locations where it is not possible, for example, due to the lack of adequate electrical power 
during eclipse periods, or due to unanticipated failures and contingencies. 

30 [0009] The drawback of the continuous firing strategy is that, for at least some portions of 
the orbit, the thruster firing is inefficient. For example, inefficiencies may occur when the 
intermediate orbit (final phase start orbit) has an apogee altitude equal to the GEO altitude 
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(synchronous altitude) and has zero incHnation. In this case, there is no useful benefit for the 
thrust provided in the region around orbit perigee, since thrusting at perigee is most useful for- 
changing the apogee altitude or inclination, and neither of these corrections are required in 
this case. 

5 [0010] To provide a higher efficiency orbit transfer, alternative strategies are possible 
where the thrusters are fired for only a portion of each orbit, in regions where firing is most 
efficient. This approach is referred to as an on/off firing strategy. For example, it is well 
known that, in the above example of an orbit with a synchronous altitude apogee, the 
thrusters may be fired efficiently for some region about apogee to correct inclination and 
10 raise the perigee altitude. One drawback of this approach, however, is that although the 

efficiency of firing can be increased, the firing can be sufficiently intermittent so that it is not 
possible to deliver the spacecraft mass to GEO orbit within the allocated time. 

[001 1] Furthermore, current on/off firing strategies tend to be "'ad hoc"". That is, the off 
and on firing regions are selected based on experience with impulsive maneuvers, such as 

15 LAE maneuvers, and preferred orbital regions for executing such maneuvers. Additionally, 
the fimctional form of the expression used to model the thrust vector trajectory also may be 
determined based on experience. The trajectory itself may then be determined using a 
numerical optimization procedure in which the trajectory model parameters are varied to 
achieve the best results. The use of an assumed thrust vector trajectory and ad hoc selected 

20 orbital firing regions imposes artificial constraints that limit the performance of these 

methods compared to methods where the optimal trajectory and on/off firing regions are 
solved for directly. 

[0012] Thus, an improved method is desired that may be used to determine thrust trajectory 
and on/off firing regions in a way that provides the improved fuel use reduction with 
25 minimum impact to the total orbit transfer time. Preferably, the method does not require the 
use of assumptions regarding the fimctional form of the thrust trajectory or the orbital 
locations where firing should occur. Also, a system is desired for implementing such a 
trajectory and firing regions on-board a spacecraft to perform a transfer from an intermediate 
orbit to a final mission orbit. 
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BRIEF SUMMARY OF THE INVENTION 
[0013] The present invention relates generally to systems and methods for transferring a 
spacecraft from a first orbit to a second orbit. In accordance with one embodiment of the 
invention, the method comprises calculating thruster-off regions within an orbit transfer in 
5 which it is efficient to turn-off spacecraft thrusters, and in those thruster-off regions, turning 
off the spacecraft thrusters. 

[0014] In accordance with another embodiment of the present invention, a method for 
efficiently transferring a spacecraft to a desired orbit comprises first computing a continuous- 
firing thrust trajectory to achieve an orbit transfer. The method ftirther comprises computing 

10 thrust effectiveness values for time intervals over the continuous-firing thrust trajectory, 

comparing the thrust effectiveness values with a thrust effectiveness threshold value, and then 
computing an intermittent-firing thrust trajectory to achieve the orbit transfer. In accordance 
with this embodiment of the invention, the intermittent-firing thrust trajectory comprises 
thruster-on regions where the thrust effectiveness value is about or above the thrust 

15 effectiveness threshold value, and thruster-off regions where the thrust effectiveness value is 
below the thrust effectiveness threshold value. 

[0015] In accordance with yet another embodiment of the present invention, the step of 
computing the intermittent-firing thrust trajectory may comprise the steps of: (1) determining 
one or more thruster-off regions for a first orbit revolution; (2) computing a first updated 

20 thrust trajectory for the entire orbit transfer using the thruster-off regions identified for the 

first orbit revolution in the calculation; (3) determining one or more thruster-off regions for a 
second orbit revolution using the first updated trajectory; (4) computing a second updated 
thrust trajectory for the entire orbit transfer using the thruster-off regions identified for the 
first and the second orbit revolutions in the calculation; and (5) continue computing thruster- 

25 off regions for each successive orbit revolution and further updated thrust trajectories until a 
final intermittent-firing thrust trajectory is determined for all orbits of the entire orbit transfer. 

[0016] In accordance with one embodiment of the invention, the thruster-on regions, the 
thruster-off regions and the final intermittent-firing thrust trajectory can be determined prior 
to carrying out the orbit transfer. 

30 [0017] In accordance with yet another embodiment of the present invention, the thrust 
effectiveness value is calculated according to the equation: 
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[0018] r(o = i-^ 

K Z 

[0019] In accordance with yet anther embodiment of the present invention, prior to the step 
of comparing the thrust effectiveness value with a thrust effectiveness threshold value, the 
method further comprises determining the thrust effectiveness threshold value. For example, 
5 the thrust effectiveness threshold value may be calculated as a function of thruster shut-off 
time, fuel savings and increase in orbit transfer time. In one embodiment, the thrust 
effectiveness threshold value is denoted and can be solved for by evaluating the integrals 

T 

T,(r,)=\y\Tdt ifr<ro 

[0020] ° where, c\ f t- i- 

^i. r| = 0 if r > 

0 

[0021]' for values of between 0 and 1 with a reasonable resolution, wherein Jj gives a 
10 relationship between the thrust effectiveness threshold value and a total increase in the 
orbit transfer time, and wherein T2 gives a relationship between the thrust effectiveness 
threshold value and a reduction in firing time. 

[0022] Yet another embodiment of the present invention comprises an orbit transfer system 
comprising thrusters and at least one controller adapted to implement the orbit transfer 
15 method as set forth herein. Still another embodiment of the present invention comprises a 
spacecraft that includes the orbit transfer systems and methods described herein. 

[0023] A more complete understanding of the present invention can be derived by referring 
to the detailed description of preferred embodiments and claims when considered in 
connection with the figures. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0024] In the Figures, similar components and/or features may have the s£ime reference 
label. Further, various components of the same type may be distinguished by following the 
reference label with a second label that distinguishes among the similar components. If only 
25 the first reference label is used in the specification, the description is applicable to any one of 
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the similar components having the same first reference label irrespective of the second 
reference label. 

[0025] Fig. 1 is a drawing of one embodiment of a spacecraft that may include the systems 
and methods of the present invention; 

5 [0026] Fig. 2 is a flow chart illustrating one embodiment of a method according to the 
present invention; 

[0027] Fig, 3 is a block diagram illustrating one embodiment of an orbit transfer system 
according to the present invention; 

[0028] Fig, 4 is graph illustrating the relationship between fiiel savings and the increase in 
10 transfer time for one embodiment of a method according to the present invention; and 

[0029] Figs. 5a and 5b are graphs illustrating the relationship between Gamma and time for 
an orbit transfer conducted with thmsters firing intermittently according to one embodiment 
of a method according to the present invention. 



1 5 DETAILED DESCRIPTION OF THE INVENTION 

[0030] The present invention relates to a orbit transfer techniques, and more particularly to 
systems and methods for efficiently transferring a spacecraft from one orbit to another orbit. 

[0031] Prior-art systems use a continuous-firing orbit transfer strategy, in which the ion 
thrusters are fired almost continuously over some period to provide a transfer from an initial 

20 to a final orbit. One drawback of firing ion thrusters continuously to carry-out an orbit 

transfer is that there may be some portions of the orbit where firing the thrusters is of only 
marginal benefit. Thus during those portions of the orbit, fuel is expended, but only a very 
small useful change in the orbit is achieved. To address this problem, the systems and 
methods of the present invention determine a spacecraft thrust vector trajectory and off-firing 

25 regions that eliminate these inefficient thrusting regions. As a result, the spacecraft orbit 
transfer requires less fuel, and thus, larger, more capable payloads can be put into space, 
increasing the revenue generating capability and utility of the spacecraft. Furthermore, the 
systems and methods of the present invention eliminate or mostly eliminate the inefficient 
firing regions, so that the increase in the orbit transfer time is minimized. 
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[0032] Referring now to Fig. 1, one embodiment of a spacecraft 100 that may include the 
systems and methods for efficiently transferring a spacecraft from one orbit to another orbit 
in accordance with the present invention is shown. In the illustrated embodiment, spacecraft 
100 includes a spacecraft center body 102, and solar panels 104a and 104b and 
5 communication antennas 106a and 106b connected to spacecraft center body 102. In 
addition, spacecraft 100 includes an attitude control system that may comprise attitude 
sensors, control torquers, and a processing system. For example, in the embodiment 
illustrated in Fig. 1, the attitude control system includes an earth sensor 108, a sun sensor 
1 10, an inertial measurement unit (e.g., a gyro) 112, reaction wheels 114, attitude control 
10 thrusters 116, and orbit transfer thrusters 118. 

[0033] As one skilled in the art will appreciate, the sensors (earth sensor 108, sun sensor 
110 and IMU 1 12) are adapted to monitor the attitude of the spacecraft, and the control 
troquers (reaction wheels 114 and thrusters 116, 118) are adapted to provide attitude control 
and orbit-change capability. This may include transferring the spacecraft from its insertion 
15 orbit to its final orbit, and maintaining the spacecraft attitude while in its final orbit. Also, as 
one skilled in the art will appreciate, the attitude control processing system can be adapted to 
control the spacecraft as it transfers from its insertion orbit to its final orbit, as well as 
controlling the spacecraft attitude in its final orbit. 

[0034] The location and position of thrusters 116 and 1 18 on spacecraft 100 are for 
20 illustration purposes only. One skilled in the art will appreciate that the thrusters can be 
located at number of different locations on the spacecraft, which can be dictated by the 
spacecraft mission and orbit transfer requirements. Also, the spacecraft illustrated in Fig. 1 is 
but one example of a spacecraft that may embody the systems and methods of the present 
invention. Therefore, the present invention is not limited to the illustrated embodiment. One 
25 skilled in the art will appreciate that the systems and methods of the present invention can be 
used on any spacecraft configuration and for any spacecraft mission. 

[0035] In accordance with the present invention, the term Gamma ( T ) will be used herein 
to refer to a thrusting effectiveness for a particular period during a given orbit transfer 
trajectory. As used herein. Gamma is the ratio of the change in the total orbit transfer time to 
30 the thruster shut off time, and may be computed for all time steps within the orbit transfer 
trajectory, as described below. At the portions of the orbit where Gamma is less than 1, it is 
beneficial to turn the thrusters off. In this case the total orbit transfer time increases, but 
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because the shut-off time is greater than the increase in orbit transfer time, the total firing 
time and fuel usage are reduced. In fact, when the thruster is shut-off at a given time step, the 
reduction in firing time, which is proportional to the fiiel savings, is just one minus the value 
of Gamma. Thus, in accordance with the present invention, the smaller the Gamma value, 
5 the greater the fiael benefit, and the lower the transfer time increase. In the optimal case 
where Gamma is equal to zero (Gamma cannot be less than zero), the transfer time does not 
increase at all when thruster firing time is removed. 

[0036] The present invention provides for systems and methods for improving the 
efficiency of trcuisferring a spacecrafl: fi-om an initial orbit to a final orbit. In accordance with 

10 the present invention, an optimal orbit transfer trajectory is solved for starting fi"om the prior- 
art continuous-firing solution. Next, the system determines when thruster can and/or should 
be tumed-off by calculating when Gamma is less than a specified threshold value. The 
threshold value is determined based on a number of factors, including desired fuel savings 
and maximum orbit transfer time available. A numerical iteration is performed to determine 

15 a minimum-time orbit transfer trajectory that achieves the final orbit with the Gamma 
constraint imposed. 

[0037] Based on the results of the numerical optimization calculations, the spacecraft 
attitude is controlled to orient the thrusters' thrust vector in the direction of the specified 
thmst trajectory vector. The thrusters then are modulated so they fire within the specified 

20 firing regions of each orbit and do not fire in the specified off regions. To carry out an orbit 
transfer, the present invention can determine parameters for thruster firings to achieve the 
desired orbital change most efficiently. The thruster parameters can include when to tum 
thrusters on and off, and other parameters which are either the thrust vector direction (in an 
inertial or any other coordinate fi*ame) or some other parameters which enable the spacecraft 

25 on-board software to compute the thrust vector direction at any time. A more detailed 
discussion of the systems and methods of the present invention is give below. 

[0038] Referring now to Fig. 2, flowchart 200 illustrates one embodiment of a method for 
carrying out orbit transfers in accordance with the present invention. In accordance with the 
illustrated embodiment, first a continuous-thruster- firing orbit transfer trajectory is calculated 
30 (block 202) using start orbit parameters 204 and mission orbit parameters 206. This 

"minimum-time" orbit transfer trajectory is first determined as if the change is to be carried 
out with continuous-firing thrusters. After the initial orbit transfer trajectory is determined, 
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the thruster effectiveness value Gamma (F ) is calculated for the initial continuous-firing 
orbit transfer trajectory (block 208). 

[0039] In accordance with one embodiment of the present invention, the thruster 
effectiveness value Gamma can be computed as a function of time during the orbit transfer 
5 according to the following Equation (1): 

[0040] r(0 = i-^ (1) 

[0041] where X is a 6x 1 vector [aJ' = [A, , , A3 , , , ]), which is the co-state vector of 
the minimum-time orbit transfer solution. The co-state vector is a function of time and varies 
at each time step over the orbit transfer trajectory. The vector X. is a Lagrange multiplier that 
10 multiplies with the constraint equations in the formulation of the minimum-time orbit transfer 
problem. As is known to those of skill in the art, the constraint equations are the equations of 
motion of the orbit elements, which are also referred to as the orbit perturbation equations. 

[0042] The orbit element state vector z is a 6x1 vector that contains equations for: (1) the 
orbit semi-major axis; (2) the y-component of eccentricity vector; (3) the x-component of 
15 eccentricity vector; (4) the y-component of inclination vector; (5) the x-component of 

inclination vector; and (6) the eccentric longitude (x and y are with respect to the equinoctial 
coordinate frame). The orbit state vector is as follows: 

a 
h 
k 
P 

q 

F 

[0044] where the orbit element perturbation equations for a,h,k, p, q, F are well known in 
20 the art. Each element of A multiplies with a corresponding orbit element perturbation 

equation. In particular, X,^ is the 6^^ element of X and multiplies the perturbation equation 
for the eccentric longitude F , which is the 6**^ element of the orbit element state vector z . In 
Equation (1), z is the 6x1 vector of orbit element time derivatives. Also, the quantity F in 
Equation (1) denotes the time rate of change of the eccentric longitude F and is also the 6* 
25 element of z . 



[0043] z = 
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[0045] After the thruster effectiveness value Gamma r(/) is calculated for the initial thrust 
vector trajectory (block 208), an on/off threshold value is determined (block 210). A 
number of factors are used to generate the value for Fq , including thruster shut-off time, fuel 
savings, and the increase in orbit transfer time, as set forth below. Once the threshold value 
5 Tq is determined, regions within revolution 1 of the orbit transfer are determined in which 
r < Tq (block 212), which determines when the thrusters will be turned off Finally, a 

minimum-time orbit transfer thrust trajectory vector is determined using the regions where 
the thrusters are turned off (block 214). 

[0046] In accordance with one embodiment of the present invention, to determine those 
10 regions in which the thrusters are to be turned off, the thruster effectiveness value Gamma 

r(/) for the initial orbit transfer trajectory is evaluated to establish a relationship between the 
Gamma threshold value , the thruster shut-off time, the fuel savings, and the increase in 
the total orbit transfer time. Once the Gamma threshold value Fq is estabUshed, it is used in 

the final optimization step that solves for the minimum-time, on/off thruster-firing orbit 
1 5 transfer thrust traj ectory vector (block 214). 

[0047] To determine the effect of the Gamma threshold value Tq , the integrals of Equation 
(2) are evaluated for all time steps along the continuous-firing orbit transfer trajectory. 

T,(r,) = ]nrdt if r < 

[004«1 ; where, ^ = Q if T > 

0 

[0049] These two integrals can be evaluated for all possible values of the Gamma threshold 
20 To varying from 0 to 1 with a reasonable resolution. The first integral T, gives the 

relationship between the Gamma threshold value and the total increase in the orbit transfer 
time. The second integral T2 gives the relationship between the Gamma threshold value and 
the reduction in firing time. The fiiel savings Mj. is proportional to the reduction in firing 
time, and in one embodiment, can be computed using Equation (3) 

T NF 

25 [0050] = ^i^-^ (3) 
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[0051] where N is the number of thrusters firing, /\ is the thrust for a single thruster, and 
is the mass specific impulse of the thruster. Using Equations (3) and (2) it is possible to 

determine the relationship between Gamma and the fuel savings. Additionally, since the 
relationship between Gamma and the increase in transfer time can be determined from 
5 Equation (2), the relationship between fuel savings and increase in transfer time can also be 
determined. Based on the relationship between fuel savings and the increase in orbit transfer 
time, and given the time available for orbit transfer, a Gamma threshold can be specified 
below which the thrusters are turned off The minimum-time orbit transfer trajectory vector 
then is computed with the constraint that the thrusters are tumed off wherever T < Tq (block 
10 214). 

[0052] The minimimi-time orbit transfer trajectory vector can be calculated using an 
iterative operation to determine the initial co-state vector A , and the minimum transfer time 
fi-om the initial orbit to the final orbit, while imposing the thruster on/off constraint of 
r < Fq . hi accordance with one embodiment of the invention, the initial co-state vector A 
15 and the minimum transfer time can be obtained using a Nev^on-Raphson iteration approach, 
in which the initial values for the co-state vector and the minimum transfer time can be the 
values which form the solution for the continuous-fire orbit transfer problem, 

[0053] After the minimum-time orbit transfer thrust vector is determined for the initial 
revolution, the regions of the next orbit revolution (revolution 2) are identified where F < Tq 

20 (block 216). The regions are determine in the manner as discussed above, but instead of 

starting the iteration firom the continuous-firing orbit transfer solution, the iteration is started 
firom the newly calculated solution with the gamma constraint applied to the first orbital 
revolution. By shutting off the thrusters in the region of revolution 2 where gamma is below 
threshold, as well as in the identified revolution 1 regions where gamma is below threshold, a 

25 new minimum-time orbit transfer trajectory is computed (block 2 1 8). This process is 

repeated for each subsequent orbit revolution until the optimal trajectory is obtained with the 
thruster off in all regions where F < Fq . 

[0054] Li accordance with an alternative embodiment of the invention, an alternative 
method is possible where the off-firing regions are first identified for all orbit revolutions of 
30 the orbit transfer. Next, the optimization problem is solved with thmster shut-off in all those 
regions. Once the solution is obtained. Gamma can be re-computed, and the thruster shut-off 
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regions may be adjusted again according to the new values of Gamma. This way the solution 
may be iterated until there is insignificant change in the solution. Regardless of the method 
used to determine the optimal orbit transfer trajectory, the solution will be such that the orbit 
transfer time is minimized, and the thruster does not fire wherever Gamma (as determined 
5 from Equation (1)) is less than or equal to the specified threshold. 

[0055] The present invention also includes a system for high efficiency spacecraft orbit 
transfer. In accordance with one embodiment of the invention, the orbit transfer system 
includes one or more processors for processing data provided from the ground, as well as data 
generated by on-board attitude sensors to carry out the method as described above. The 

10 processor(s) are operatively connected to the attitude sensor(s) and to command receivers and 
antennas to receive the data. To provide attitude control and effect an orbit change, the 
processor(s) are operatively connected to thrusters, thruster gimbals, solar array drives, and 
/or reaction wheel assemblies. As described above, an inertial thrust vector trajectory for the 
orbit transfer is generated. The system is adapted to orient the spacecraft so that the 

15 spacecraft thruster thrust vector remains aligned with the inertial thrust vector trajectory. The 
spacecraft thrusters are turned on and off during portions of each orbit until the spacecraft 
transfers to the desired orbit. 

[0056] Referring now to Fig. 3, one embodiment of an orbit transfer and attitude control 
system 300 in accordance with the present invention is shown. System 300 comprises 
20 attitude determination logic 302, an ephemeris generator 304, a trajectory and thruster on/off 
command generator 306, an attitude command generator 308, a maneuver supervisor 3 1 0, 
attitude and solar array control logic 312, and actuator command logic 314. 

[0057] In accordance with one embodiment of the present invention, attitude determination 
logic 302 is adapted to maintain an estimate of the spacecraft inertial attitude by propagating 
25 three-axis gyro data from an Inertial Measurement Unit (IMU). The inertial attitude and gyro 
bias estimates can be updated utilizing attitude sensor data, including earth and sun sensor 
data or star tracker data. 

[0058] Ephemeris generator 304 is adapted to determine the spacecraft ephemeris from 
ground supplied parameters, and other information. To generate the ephemeris, ephemeris 
30 generator 304 can be adapted to propagate the orbit perturbation equations using estimates of 
the applied thruster force and perturbation forces due to the sun and/or moon, and other 
higher order earth gravitation effects. 
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[0059] Trajectory and thruster on/off command generator 306 is adapted to generate the 
thrust vector trajectory using ground supplied parameters and information from the on-board 
ephemeris generator. Additionally, trajectory and thruster on/off command generator 306 can 
be adapted to generate thruster firing information that indicates when the thrusters should be 
5 tumed on and when they should be turned off. Both the thrust vector trajectory and firing 
state may be generated based on parameters uploaded to the spacecraft. The ground supplied 
parameters can include parameters representative of the on/off thruster firing orbit transfer 
trajectory determined using the invented method, as well as other information. 

[0060] In accordance with one embodiment of the present invention, ephemeris generator 
10 304 is adapted to integrate both the orbit perturbation equations and the co-state differential 
equations to obtain both the ephemeris parameters and the co-states at a given time. 
Trajectory and thruster on/off command generator 306 then can use the co-states to solve for 
the thrust vector direction. The orbit perturbation equations, the co-state equations, and the 
equation relating the co-states to the thrust vector direction are the same equations used to 
15 generate a continuous-firing minimum-time orbit transfer solution, and hence are well known 
to those skilled in the art. 

[0061] Trajectory and thruster on/off command generator 306 also can generate thruster 
activation information based on ground-supplied parameters and the spacecraft ephemeris. 
As discussed above, thruster on/off commands can be determined by using the calculated 
20 Gamma value, which is based on the ephemeris and co-state information according to 

Equation 1. When Gamma is above the specified Gamma threshold value, the thrusters are 
fired, and when Gamma is below the specified Gamma threshold value, the thrusters are 
tumed off 

[0062] Trajectory generator 306 calculates an attitude trajectory using the generated thrust 
25 vector trajectory. This attitude trajectory can be represented as a coordinate transformation 
from an inertial to a spacecraft body coordinate frame. When the spacecraft tracks this 
attitude trajectory, the thruster thrust vector remains aligned with the thrust trajectory vector. 
In one embodiment of the invention, the thrust vector may be close to the yaw(x)/roll(y) 
plane of the spacecraft body, for example, close to the roll axis of the spacecraft. 
30 Furthermore, the attitude trajectory may be determined so that the spacecraft is controlled to 
an orientation about the thruster thrust vector. In this orientation, the sun vector may Ue 
almost in the spacecraft yaw/roll plane. Since the solar array rotates about the pitch(z) axis 
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(the normal vectors to the solar array panel surfaces are in the yaw/roll plane), the solar 
arrays may be oriented to point at the sun while the orbit transfer is performed. The solar 
arrays then can provide adequate power to the spacecraft to operate the ion thrusters and 
other on-board spacecraft equipment. 

5 [0063] The thruster firing information from trajectory and thruster on/off command 
generator 306 is supplied to maneuver supervisor 310, which controls the initiation and 
completion of orbit transfer thruster firing segments. Maneuver supervisor 310 can switch 
operating modes, switch actuators, and/or determine which attitude the spacecraft tracks 
during periods when the thrusters are not firing. 

10 [0064] The spacecraft attitude command generator 308 receives inertial attitude 

information fi-om attitude determination logic 302 and attitude trajectory information fi"om 
trajectory generator 306 and generates attitude error correction information. The attitude 
error information is computed based on the spacecraft inertial attitude, the thrust vector 
trajectory, and possibly the Sun vector, and is used to null the attitude errors so that the 

15 thruster thrust vector remains aligned with the thrust vector trajectory, and the sun vector 
remains oriented in the spacecraft yaw/roll plane (perpendicular to the solar array rotation 
axis). 

[0065] The attitude error information fi-om attitude command generator 308 is input to 
attitude and solar array control logic 312, which generates solar array position commands 

20 necessary to point the solar array at the sxm. While the thrusters are firing, attitude and solar 
array control logic 312 can partition the torque demand between the thruster gimbals and the 
reaction wheel assemblies (RWAs). When the thrusters are not firing, the attitude and solar 
array control logic 312 applies the control torque demand entirely to the RWAs, and/or 
partitions the control torque demand between the RWAs and hydrazine control thrusters 

25 (REAs). 

[0066] In response to the commands generated by attitude and solar array control logic 312, 
actuator command logic 314 sends commands to the RWAs, the thruster gimbals, the ion 
thruster gimbals, and/or the solar array drives, which in turn controls the attitude of the 
spacecraft. Actuator command logic 314 also may control the RWAs, thruster gimbals, ion 
30 thruster gimbals and/or solar array drives using information and commands fi-om maneuver 
supervisor 310. 
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[0067] Referring now to Figs. 4-6, an example of the benefits of the systems and methods 
of the present invention will be described. This particular example assumes an orbit transfer 
using two Hall Current Thrusters (HCTs) fired at about 9 kW. The starting orbit has perigee 
and apogee altitudes of about 13,098 km and about 35,786 km, respectively, and 7 degrees 
5 inclination. The initial spacecraft mass is about 4500 kg, and the desired final orbit is a 

circular Geosynchronous orbit with zero inclination. Fig. 4 is a graph 400 illustrating the plot 
between the firing time reduction, which is proportional to the fiiel savings (y-axis), and the 
increase in total transfer time (x-axis). From Fig. 4, as an example, it can be noted that by 
allowing an additional 5 days to transfer the spacecraft to its final orbit, about 5.5 days of 
10 firing time can be removed. This point on the curve shown in Fig.4 corresponds to a gamma 
threshold value of 0.6, and during the orbit transfer the thrusters are shut off wherever the 
value of Gamma is less than or equal to 0.6. The 5.5 days of reduced thmster firing 
corresponds to a fiiel savings of roughly 15.7 kg, assuming a total HCT thrust of 0.13 lbs 
(with two HCTs firing) and an HCT specific impulse of 1771 seconds, 

15 [0068] Fig 5a shows a graph 500 of Gamma versus time during an entire orbit transfer 
using the prior-art continuous-firing approach. Fig 5b is an expanded version of Fig 5a for 
the first 10 days of the orbit transfer. The transfer time using the continuous-firing approach 
is 81 .1 days, which is the minimum possible time for this orbit transfer using two HCTs and 
the HCT, spacecraft, and orbit parameters given previously. 

20 [0069] As per the approach shown in Figure 2, the shut-off region in the first revolution is 
identified based on Gamma less than 0.6 from Figure 5b. Next, the minimum time 
continuous solution with the first revolution off constraint is obtained. A similar Gamma 
versus time relation plot for this new solution is obtained. Next, the shut-off region in the 
second revolution is identified where Gamma is less than 0.6 in the new plot. This process is 

25 then repeated for 37 revolutions. The total transfer time increase obtained is 4.4 days and the 
firing time reduction is 3.9 days, which corresponds to a fiiel reduction of 1 1.1 kg. 

[0070] In conclusion, the present invention provides novel systems and methods for 
efficiently transferring a spacecraft from an initial orbit to a final orbit. While detailed 
descriptions of one or more embodiments of the invention have been given above, various 
30 altematives, modifications, and equivalents will be apparent to those skilled in the art without 
varying from the spirit of the invention. For example, while various embodiments for solving 
for the orbit transfer thrust vector trajectory and thruster on/off times are given herein, one 
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which is defined by the appended claims. 
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